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1 Introduction 



The study of primordial chemistry of molecules adresses a number of interesting 
questions pertaining to the thermal balance of collapsing molecular protoclouds. 
In numerous astrophysical cases molecular cooling and heating influence dynam- 
ical evolution of the medium 

The existence of a significant abundance of molecules can be crucial on the dy- 
namical evolution of collapsing objects. Because the cloud temperature increases 
with contraction, a cooling mechanism can be important for the structure forma- 
tion, by lowering pressure opposing gravity, i.e. by allowing continued collapse 
of Jeans unstable protoclouds. This is particularly true for the first generation 
of objects. 

It has been suggested that a finite amount of molecules such as H 2 , HD and LiH 
can be formed immediately after the recombination of cosmological hydrogen 
(Lepp and Shull 1984, Puy et al 1993). The cosmological fractional abundance 
HD/H2 is small. Nevertheless, the presence of a non-zero permanent electric 
dipole moment makes HD a potentially more important coolant than H 2 at 
modest temperatures, although HD is much less abundant than H 2 . Recently 
Puy and Signore (1997) showed that during the early stages of gravitational col- 
lapse, for some collapsing masses, HD molecules were the main cooling agent 
when the collapsing protostructure had a temperature of about 200 Kelvins. 
Thus in section 2, we analytically estimate the molecular cooling during the 
collapse of protoclouds Then, in section 3, the possibility of thermal instability 
during the early phase of gravitational collapse is discussed. 



2 Molecular functions of a collapsing protocloud 

We consider here only the transition between the ground state and the first state. 
The transition is a dipolar transition whereas for H 2 molecule the transition is 
quadrupolar (because it has no permanent dipolar moment). 
In this study, we consider only the molecular cooling which is the dominant 
term in the cooling function as we have shown for a collapsing cloud (Puy & 
Signore 1996, 1997, 1998a, 1998b). we obtain for the cooling function of HD 
the expression: 
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and the molecular cooling for the H 2 molecule is given by: 
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The total cooling function is given by 

A-Total = Afffl + Aff 2 = Ahd(1 + £h 2 ) 
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where 
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is the ratio between the H2 and HD cooling functions. Finally we deduce: 
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and for the ratio: 
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We study a homologous model of spherical collapse of mass M similar to 
the model adopted in Lahav (1986) and Puy & Signore (1996) in which we only 
consider H2 and HD molecules. The aim of this paper is to analyse, through 
our approximations (where the only two first excited levels are considered), the 
evolution of the cooling and the potentiality of thermal instability. 
First, let us recall the equations governing dynamics of a collapsing protocloud: 

dT m _T dr 2 

^Total (7) 



T dr 
-2-.— 

r dt 



(8) 



dt r dt 3nk 

for the evolution of the matter temperature T m , 

d 2 r _ 5kT m GM 
dt 2 1mnr r 2 
for the evolution of the radius r of the collapsing cloud and 

dn dr 
dt r dt 

for the evolution of the matter density n. 

At the beginning of the gravitational collapse the matter temperature increases, 
then due to the very important efficiency of the molecular cooling, the temper- 
ature decreases. We consider here this transition regime i.e. the point where 
the temperature curve has a horizontal asymptot (see Puy & Signore 1997): 
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Finally we obtain for the total molecular cooling 
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M 2 with k = 1.4 X 10 22 , T = 128.6 K, AT = 383.4 K and 
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3 Thermal Instability 



We have learned much about thermal instability in the last 30 years since the 
appearance of the work of Field (1965). Many studies have focused on the 
problem of thermal instabilitie in different situations. In general, astronomical 
objects are formed by self-gravitation. However, some objects can not explained 
by this process. For these objects the gravitational energy is smaller than the 
internal energy. Thus if the thermal equilibrium of the medium is a balance 
between energy gains and radiative losses, instability results if, near equilibrium, 
the losses increase with decreasing tempraturc. Then, a cooler than average 
region cools more effectively than its surroundings, and its temperature rapidly 
drops below the initial equilibrium value. 

Thus, if we introduce a perturbation of density and temperature such that 
some thermodynamics variable (pressure, temperature...) is held constant, the 
entropy of the material S will change by an amount 5S, and the heat-loss 
function, by an amount 5^f = T — A. From the equation 



there is instability only if: 
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We know that TdS = C p dT in an isobaric perturbation. The corresponding 
inequality can be written: 

\St) p = \5f)n + df\5^)T < 

which characterizes the Field's criterion. The criterion involves constant P be- 
cause small blobs tend to maintain pressure equilibrium with their surroundings 
when heating and cooling get out of balance. In our case (i.e. at the transition 
regime), a thermal instability could spontaneously be developped if the Field's 
Criterion is verified: 

dink 
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In our case the logaritmic differential of the total molecular cooling function is 
given by 
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In figure 1, we have plotted the curve y(T m ) = dink/ dlnT m . It shows that 
the Field's criterion is always verified. 

We conclude that a thermal instability is possible at the transition regime. Thus 
by maintcning the same pressure in its surroundings, such a blob would get 
cooler and cooler and denser and denser until it could separate into miniblobs. 
This possibility is very interesting and could give a scenario of formation of 
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Figure 1: dlnATotai/dlnTm in y-axis and temperature of the matter in x-axis (in 
Kelvins) 



primordial clouds. However, a quantitative study is necessary to evaluate the 
order of magnitude of the mass and the size of the clouds. This last point is 
crucial. 

Therefore, from these estimations, the calculations of molecular functions could, 
in principle, be extended to the CO molecules for collapsing protoclouds at 
z < 5 (see Puy & Signore 1998b). But, in any case -and in particular for the 
discussion of thermal instability at the transition regime- the corresponding total 
cooling function must be numerically calculated because the excitation of many 
rotational levels must be taken into account. This extended study is beyond the 
scope of our paper. 
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